Formation of energy gap in higher dimensional spin-orbital liquids 
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A Schwinger boson mean field theory is developed for spin liquids in a symmetric spin-orbital 
model in higher dimensions. Spin, orbital and coupled spin-orbital operators are treated equally. 
We evaluate the dynamic correlation functions and collective excitations spectra. As the collective 
excitations have a finite energy gap, we conclude that the ground state is a spin-orbital liquid with 
a two-fold degeneracy, which breaks the discrete spin-orbital symmetry. Possible relevence of this 
spin liquid state to several realistic systems, such as CaV4Vg and Na2Sb2Ti20, are discussed. 
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The formation of a spin gap in two- or higher dimen- 
sional quantum spin systems is a long-standing issue in 
strongly correlated problems. jO Several physical mech- 
anisms were proposed to explain the formation of the 
spin gap in the low-energy excitations. Most of them 
focus on one-dimensional spin chains and spin ladders, 
such as an S=l antiferromagnetic Heisenberg chain and 
even-leg spin ladders. In higher dimensions Anderson 
proposed that strong quantum fluctuations for spin-1/2 
systems may destroy the antiferromagnetic long-range 
order in two dimension, and lead to form a resonat- 
ing valence bond (RVB) state where a spin gap may 
open. |Q However, it becomes true only for some frus- 
trated spin systems such as on the Kagome lattice or in 
the Majumdar-Ghosh model with a strong next nearest 
neighbor interaction, otherwise there exists antiferromag- 
netic long-range orders in the ground state on a square 
and triangle lattice. Recently it has been realized that 
orbital degrees of freedom of d- and f-electrons in transi- 
tion metal ions provide a new route to find the physical 
mechanism of spin gap formation. Several spin-orbital 
models have shown the tendency of the formation 

of a spin gap in the ground states due to strong orbital 
and spin quantum fluctuations. Accumulating numerical 
calculations show that spin liquid state may be formed in 
some one-dimensional spin-orbital coupled systems. Be- 
haviors in higher dimensional systems are relatively less 
clear. There are several higher dimensional spin-gap ma- 
terials such as Na2Ti2Sb20 || and NaV40g flfj] , in which 
the orbital degree of freedom might play a key role in the 
formation of the spin gap. 

The spin-orbital model Hamiltonian is written as 

HI 
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(i) 



This model may be derived from an electronic model 
with double orbital degeneracy. This model possesses 
an SU{2) <g> SU{2) symmetry. When J s = J T = J , the 



model have an additional discrete symmetry, namely, the 
permutation symmetry between spin and orbital. In this 
case the model can be written as a combination of two 
symmetric models 



H = (J + Jo) g [2Si ■ Si + I) [2n . Tj + \ 

-J J2(2S i -S j -l)(2n-T j -l). (2) 



The former part is the standard SU(4) spin-orbital 
model, which is solvable in one-dimension and has been 
investigated intensively. || The second part is the model 
first proposed by Santoro ||, and also possesses the 
SU(4) symmetry with different generators. However the 
combination of the two models breaks the SU(4) symme- 
try. 

^From previous studies, it has been found that the in- 
terplay between spin and orbital degrees of freedom pro- 
duces either quantum ordered or disordered phases. Spin 
liquid states with an energy gap were found in one dimen- 
sions. To investigate the model systematically, we 
try to develop a simple theory, which can describe the dis- 
ordered state with an energy gap as well as the ordered 
states. The Schwinger boson theory is an ideal candidate. 
The theory was first used to the spin SU(2) Heisenb erg 
model, and then was generalized to SU(N) systems. |l3|] 
The advantage of this theory is that it can describe either 
quantum ordered or disordered states. The results for 
two- and three-dimensional Heisenberg model are con- 
sistent with the spin wave theory very well. Here we 
present the Schwinger boson mean field theory for this 
spin-orbital system. 

For the present model, there are four possible states 
on each site i according to the eigenvalues of Sf and 
7f: |1) = | + l/2, +1/2), |2) = |-l/2, +1/2), |3) = 
1 + 1/2,-1/2), |4) = |-l/2,-l/2). We introduce four 
Schwinger bosons to describe these four states: \fi) = 
a tn |0) j where |0) is the vacuum states and /i — 1,2, 3, 4. 

There is a constraint for the four bosons, Y) „— x a lp a v = 
1, for each site. On these basis the spin and orbital oper- 
ators can be expressed in terms of these four Schwinger 
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bosons. 

&t = a li«i2 + a\ 3 ai4; t+ = a^aa + a\ 2 a,iA) 
S i = a i2 a ii + a\ 4 a i3 ; rf = <4 3 <Xii + aj 4 a i2 ; 

S'i = ^ ( a li fl ii ~ a h a *2 + a l3 a i3 - a^cii^J ; 

r f = ^ ( a li a ii - aj 3 a l3 + aj 2 a i2 - al 4 a i4 

Thus, the Hamiltonian is rewritten in terms of the 
Schwinger bosons 
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where b{k) — — 2zJo(Af 4 — Af 3 )7 e (fc), The Kronecker 
product for block matrices is used, um The Hamiltonian 
can thus be diagonalized analytically, and the free energy 
is evaluated to establish the mean field equations. Due 
to the symmetry in the Hamiltonian, there exists two 



sets of solutions: (I) A 44 
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With the notation 



A, Ai, and A 2 , the same branches of spectra are given 
by 



where A ij4iv = a t ^a 0l/ - a lu a m and By )Atv = a-i^a^ + 
a-ivCijfj,. We have introduced antisymmetric and symmet- 
ric operators A and B for the purpose of the mean field 
calculations. The following theory is limited to the case 
—J > Jo > 0. We should introduce different order pa- 
rameters in the different parameter range. The local La- 
grangian multiplier is introduced to realize the local con- 
straint for hard core bosons. In the mean field approach 
we shall take it as site-independent A. The thermody- 
namic averages of the operators A and B are introduced 
as the order parameters, respectively, 

(Ai jjMV ) = -2iA° (n - rj ); {B ijjMV ) = 2A« (r< - r 3 ). 
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j) and A e ^ u (ri — rj) are odd and even functions 



with respect to the indices fj,, v or the sites r^, rj. In the 
momentum space, we take 
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where 5 points to the nearest neighbor sites. By utilizing 
the Pauli matrices a a (a = x,y,z ) and the identity 
matrix er , the decoupled Hamiltonian can be expressed 
in a compact form of 8 x 8 matrix, 



I 
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where 



*l = (al 1 ,al 2 ,al 3 ,al 4 ,a-ki,a- k 2,a-k3,a- k 4); 
H(k) = Acr ® cr (g) cr — i&y ® A{k) + b(k)a x ® cr. 
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where a{k) = —2z(J + J Q )Aj s (k) and b±(k) = -2z(J + 
J )Ai7 s (fc) ± 2zJ A 2 7 c (fc). Moreover, these two ground 
states are degenerated. The degeneracy originates from 
the symmetry of permutation of spin and orbital opera- 
tors. When Jo = 0, the model is reduced to the standard 
SU(4) spin-orbital model. The two spectra become de- 
generated, 



w(A) = yJ\*-[2z(J+ Jo) 7s (fc)] 2 (A 2 2 + A 2 3 



This way we recover the spectra for the SU(N=4) model. 

In the following we will focus on the ground state with 
an energy gap. So we do not consider the Bose-Einstein 
condensation which may give rise to the long-range or- 
der. Indeed the long-range orders appear when Jo = 
or J — — Jq on two- and three-dimensional hypercubic 
lattices. [[ll| The condensation occurs when the system 
deviates from those symmetric point slightly. General 
description of complete solutions to the problem will be 
presented elsewhere. To determine the order parame- 
ters, we introduce a set of dimensionless parameters A = 
2z(J+ J ) A/A, Ai = 2z(J+ J )Ai/A, A 2 = 2zJ A 2 /A, 
and A = X/[z( J + Jo)], then the quasiparticle excitation 
spectrum becomes 



M*) = V 1 - A2 7, 2 (fc) - [Ai7.(fc) ± A 27c (fc)] 2 . (5) 
The self-consistent mean field equations are given by 
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Substituting Eq.(^) into Eq.(||), we have 



dk 
(2n) d 
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A 2lc {k) ls {k) = 0. 



If A 2 7^ 0, Ai must be equal to zero. Oppositely, if 
A 2 = 0, the solution is for the case of Jo = 0. Therefore 
for Jo > 0, the solution is Ai = with <D+(fc) = = 
Co{k). The two spectra are also degenerated. Two sets of 
saddle point solutions become, corresponding to the spin 
liquid phase with an energy gap in elementary excita- 
tions: (I) A? 4 - Af 3 = A 2 , A° 12 = A° 34 = A, and A°„ = 
otherwise. (II) Af 4 - Af 3 = A 2 , AJ 3 = -A^ 4 = A, and 
A° = otherwise. We focus on the first set of solu- 
tions and then present the results for the second set of 
solutions. 

The one-particle Green's function in an 8 x 8 matrix 
form is defined as 

G(M) = -*(0|T{$ fc (i)$t(0)}|0 

where |0) is the ground state of the Hamiltonian. Its 
Fourier transform is given by 

G(k,u>) = ((a; + i8)a z <g> ct <S> cr - Hu)^ 1 
1 

~ uj 2 -X 2 + a 2 (k)+b 2 (k)+iS X 

uia z ® Co ® °o — Acto ® fo ® a o 
+a(k)a y cto (K> cr y — b(k)a x ® a y ® a y 

with a(&) = -2z(J + J )Aj s (k) and fe(fc) = 
— 2z JoA 2 7c(A;). The corresponding saddle point equa- 
tions then become: 
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For a given value of J/ Jo, we have a set of solutions 
for A, A, A 2 . Our solutions in this paper are limited to 
min(a)(fc)) 7^ 0. To determine the physical properties, we 
evaluate the dynamic correlation functions for the spin 
Sf, orbital Tf, and spin-orbital density operators 2SfT[. 
After some algebra, we have 



1 f dk t C x {k,g) 
' 8 J {2ir) d \oj(k)uj(k + q) 
1 1 



Cl + id + Lu(k) + uj(k + q) Cl + iS- uj(k) - uj(k + q) 



with {X = S, T, ST) 

C s (k, q) = \ 2 -a (k) a(k + q) - b(k)b(k + q): 
C T (k, q) = X 2 +a (k) a(k + q) - b(k)b(k + q): 
C ST (k, q) = X 2 -a (k) a(k + q) + b(k)b(k + q). 

There are a set of relations 
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0; 

xs T (Q,ny, 

XST(0,ty; 

XT(o,ty + XT(Q,n) 

XST(0,ty + XST(Q,n), 



1*^(0,0)] = \\ ^ 
ImIm[x T (Q,fi)] = ■ 



dk a 2 (k) 
dk b 2 {k) 



{2TT) d U 2 {k) 
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If the minimum of u>(k) is non-zero, Im[xx(q, ^)] be- 
come non-zero ONLY when 57 > 2min(w(fc)). Thus the 
collective excitations for the density-density correlation 
function have a finite energy gap, A gap = 2 min(w(fc)). 
It is worth mentioning that the solution has broken 
the discrete the permutation symmetry of spin and or- 
bital. This can be seen from the fact that, in general, 
Xs(q,£l) 7^ Xt(<Z)^)- The same expressions are obtained 
for the second set of solutions if we permute the indices S 
and T. The spectra and free energy as well as the energy 
gap are identical to the first set of solutions. Thus the 
two sets of solutions are energetically degenerated. The 
symmetries in the two states are different. More impor- 
tant, the double degeneracy of the ground state was also 
observed in one-dimensions in other approaches. [p"6p|] 
Therefore, this two-fold degeneracy is not a consequence 
of the mean field approaches and can be regarded as an 
evidence to support our mean field theory. 

Now we come to evaluate the energy gap by solving the 
mean field equations. On a one-dimensional chain, the 
energy gap can be evaluated analytically by introducing 
a parameter xq = (A 2 — A 2 )/(l — A|). The energy gap 
and the ratio of Jo / J are 



Ag a p = 4z(J+ J ) X 



TTXq 

-K(x )-E(x ) 



£ = , , K(xq)-E(x ) 
Jo E(x ) - (1 - x )K(x ) 

where K(x) is the complete elliptic integral of the second 
kind and E(x) is the complete elliptic integral of the first 
kind. We have established a one-to-one correspondence 
between the ratio Jo/ J and the energy gap. We find that 
there is a turning point at Jo/ J = 1 . The theory fails at 
the symmetric point J = 0. The energy gap still opens 
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of these materials may be relevant to the spin-orbital 
liquid with the energy gap as we discuss in this paper. 

In conclusion we develope a Schwinger boson mean 
field theory for the formation of a spin-gap in any dimen- 
sions. The ground state of the model is a spin-orbital liq- 
uid with an energy gap in an extensive parameter regime. 
This ground state breaks the discrete symmetry of per- 
mutation of spin and orbital, and is doubly degenerated. 

This work was supported by a RGC grant of Hong 
Kong and a CRCG grant of the University of Hong Kong. 



FIG. 1. The energy gap via the ratio of Jo/ J in two- and 
three- dimensions. 

The energy gaps for two- and three-dimensional lat- 
tices are plotted in Fig.l. We find that the energy gap 
opens in the regime of 0.380 < Jq/J < 9.84 for d = 2 and 
of 0.666 < Jq/J < 1.667 for d = 3. The gap closes at 
two critical ratios Jq/J. From Fig.l, it is shown that the 
energy gap appears in a larger parameter range in two- 
dimension than in three dimension.. This is consistent 
with the fact that the quantum fluctuations are stronger 
in two dimension. Out of the above parameter regimes 
the Bose condensations have to be considered, otherwise 
the mean field equations have no solutions. In this case 
the ground state may possess long-range orders as we 
discussed in the symmetric point J = — Jo- (u| We shall 
discuss these phases elsewhere. As far as we know the 
Schwinger boson mean field theory is very successful for 
the spin liquid state for s = 1 in one dimension, and anti- 
and ferromagnetic states for higher dimensions. Our the- 
ory shows it also works very well for spin-orbital liquid 
states in higher dimensions. 

The formation of a spin energy gap indicates that the 
ground state is a spin-orbital liquid. Experimentally the 
energy gap can be measured from magnetic susceptibility. 
There are several higher-dimensional materials such as 
Na 2 Ti 2 Sb 2 and CaV 4 9 . |o| It is believed that the 
orbital degrees of freedom plays an important role in the 
formation of spin gap. |0,|l8[ The low temperature phases 
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